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ABSTRACT 



Context. In spite of a large number of recent publications about the massive stellar cluster Westerlund 1, its distance from the Sun 
remains uncertain with values as low as 1.1 kpc, but largely between 4 and 5 kpc. 

Aims. The goal of this study is to determine a distance to Westerlund 1 independent of the characteristics of the stellar population and 
to study its neutral environment, using observations of atomic hydrogen. 

Methods. The H i observations are taken from the Southern Galactic Plane Survey to study H i absorption in the direction of the H n 
region created by the members of Westerlund 1 and to investigate its environment as observed in the H i line emission. A Galactic 
rotation curve was derived using the recently revised values for the Galactic centre distance of R Q = 7.6 kpc, and the velocity of the 
Sun around the Galactic centre of O = 214 km s _I . This rotation curve successfully predicts the location of the Tangent point gas and 
the velocity of the Sagittarius Arm outside the solar circle on the far side of the Galaxy to within 4 km s _I . Compared to the typically 
used values of R Q = 8.5 kpc and O = 220 kms -1 this reduces kinematically determined distances by more than 10 %. 
Results. The newly determined rotation model leads us to derive a distance of 3.9±0.7 kpc to Westerlund 1, consistent with a location 
in the Scutum-Crux Arm. Included in this estimate is a very careful investigation of possible sources of error for the Galactic rotation 
curve. We also report on small expanding H i features around the cluster with a maximum dynamic age of 600,000 years and a larger 
bubble which has a minimum dynamic age of 2.5 million years. Additionally we re-calculated the kinematic distances to nearby H n 
regions and supernova remnants based on our new Galaxic rotation curve. 

Conclusions. We propose that in the early stages of the development of Wd 1 a large interstellar bubble of diameter about 50 pc 
was created by the cluster members. This bubble has a dynamic age similar to the age of the cluster. Small expanding bubbles, with 
dynamical ages ~ 0.6 Myr are found around Wd 1, which we suggest consist of recombined material lost by cluster members through 
their winds. 
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;> 1. Introduction 



Westerlund 1 (Wd 1) is a highly reddened compact cluster with 
a large population of post-main sequence massive stars, in- 
cluding OB supergiants and hypergia nts, red and yellow super- 
giants, and Wolf- Rayet (WR) stars dClark & Negueruelall2002l: 
IClark et al.1 l2005|). The mass of Wd 1 is likely to be in ex- 
cess of 10 5 M o dClark et aUl2005l) . exceeding that of any of the 
other known massive Gal actic clusters, including NGC 6303 
dCrowther & Dessartl[T998|), t he Arches dFiger et al J 120021) and 
Quintuplet dFiger et alj|1999l) . and is more comparable to the 
mass of Super Star Clusters (SSC), previously identified only 
in other galaxies. If Wd 1 is indee d a Super Star Cluster within 
our own Galaxy dClark et alj|2005l) . this is a unique opportunity 
to study the properties of a nearby SSC, where it is possible to 
resolve the individual massive stars, and determine basic prop- 
erties more readily than in the typically more distant examples. 

In spite of a large number of recent observations of Wd 1 , its 
distan ce from the Sun remains somewhat uncertain. Piatt i et all 
(1998) derived a distance of 1.1 ± 0.4 kpc based on photome- 
try of the OB supergiants identified at the time. Using the yel- 
low hypergiants where a spectroscopic luminosity discriminant 
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is available that is lacking for OB supergiants. lClark et alJ d2005h 
determ ined an upper limit for the distance of 5.5 kpc. lClark et al.l 
(2005) also argue for a lower limit on th e distance of 2 kpc based 
on an analysis of radio continuum data dClark et alJI 19981). They 
also a rgue that the comparatively low distance of iPiatti et al.l 
(1998) is the result of an error in both the absolute magnitude 
calibration and th e reddening law that was used. Additionally, 
Clar k et aT] d2005l) did not identify any dwarf or giant stars, with 
all the OB stars being supergiants, unlike the work of Piatti et 
al. where detection of the main sequence was claimed. Further 
support for a distance in excess of 2 kpc come s from analysis 
of ne ar-IR photometry of the WR stars in Wd 1 dCrowther et al.l 
2006), where the average distance modulus derived from 23 of 
the 24 known WR stars leads to a distance of 4.7 ±1.1 kpc. The 
relatively large uncertainty in this estimate is a result of the un- 
certainty in the absolute magnitude calibration of WR stars. A 
more precise estimate of distance comes from an initial analysis 
of deep IR imaging with the VLT that reveals the main sequence 
and pre-main sequence populations in Wd 1 fro m which a pho- 
tomet ric distance of 4.0 ± 0.3 kpc is deduced (iBrandner et al.l 
120051) . 

The most recent analyses appear to be converging on a dis- 
tance in the range 4-5 kpc. In this paper, we examine H i data in 
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Fig. 1. The Galactic rotation curve toward Galactic longitude of 
339°6 using a Galactocentric radius of R Q = 7.6 kpc and a cir- 
cular velocity of = 214 km s _1 for the Sun. We assumed a 
flat rotation curve and purely circular rotation. The approximate 
locations of the spiral arms in radial velocity are indicated. Here 
the indices I and II indicate the first and second time the line of 
sight passes through that spiral arm (see also Fig.|2|. 

the direct ion of Wd 1 from the Southe rn Galactic Plane Survey 
(SGPS) (iMcClure-Griffiths et ai1l2005l) to estimate a distance to 
Wd 1 based on Hi absorption and related Hi features. We also 
examine the large-scale distribution of H i in the environment of 
Wd 1 to search for evidence of the impact of the putative SSC 
on the surrounding interstellar medium. 

2. The SGPS H i data set 

The S outhern Galactic Plane Survey (IMcClure-Griffiths et al.l 
120051) is part of the International Galactic Plane Survey (IGPS), 
a project devoted to mapping the neutral hydrogen in the plane of 
the Galaxy, the Milky Way, at arcminute scale. The SGPS maps 
the region of the Milky Way visible in the southern hemisphere. 
Other p arts of the IGPS ar e the Canadian Galactic Plane Survey 
(CGPS, Tayl or etalJ2003h . which covers the Galactic plane seen 
in the n orthern hemisph ere, and the VLA Galactic Plane Survey 
(VGPS, IStil et al.1 12006). which connects the two other surveys 
through the first quadrant of the Milky Way. 

The SGPS data were obtained with the Australia Telescope 
Compact Array (ATCA) and the Parkes 64-m radio telescope. 
The survey covers an area of 325° 2 from 253° to 358° and 5° 
to 20° in Galactic longitude and from — V.5 to + 1°5 in Galactic 
latitude. In the direction of Wd 1 the total velocity coverage is 
600 km s _1 from -300 to +300 km s _1 . The angular resolution 
is ~ 2' and the sensitivity ~ 1 .6 K. 

3. Results 

3.1. Galactic kinematics in the direction of Westerlund 1 

We carefully investigated the dynamics along the line of sight to- 
wards Wd 1 . As a starting point we used a flat rotation model to 
describe the circular rotation in our Galaxy. For this model we 




Fig. 2. The spiral arm pattern of the Milk y Way Galaxy 
(Tayl or & Cordei 1 1993b ICordes & Lazlol 12002). The names of 
the spiral arms and the line of sight (LOS) toward Wd 1 are in- 
dicated. 



have to determine the Galactocentric radius of the Sun R Q and 
its circular velocity around the Galactic centre © . Instead of as- 
suming the standard values, R Q = 8.5 kpc and = 220 kms -1 , 
which are still used i n most kinematic dista nce estimates and 
Galactic models (e.g. ICordes & L azio 2002|), we carefully in- 
vestigated the literature to find the m ost recent estima t es. i? 
was determined to be 7.6 ± 0.3 kpc by lEisenhauer et all (2005) 
from a high precision measurement of the three dimensional or- 
bit of the star S2, which is orbiting the black hole at the cen- 
tre of the Galaxy at a distance of less than 0.2". This is sig- 
nifica ntly smaller (~ 11 %) tha n the previously assumed dis- 
tance. iFeast & W hitelock (t 19971) determined the angular veloc- 
ity of circular rotation at the position of the Sun Q Q , from the 
Oort constants using data from the Hipparcos satellite. O = 
fa = 27.2 + 0.8 km s~' kpc" 1 gives = 207 ± 10 km s _1 . 
Another technique t o determine Q Q is the p roper motion of 
the Galactic Centre. iReid & Brunthalerl d2004l) found a proper 
motion of 6.4masyr~' for Sagittarius A, which translates to 
a value of 229 ± 10 km s _1 for a distance of 7.6 kpc to the 
Galactic C entre. Correcting th i s for t he peculiar motion of the 
Sun, which IReid & Brunthalerl (|2004|) suggest is between 5 and 
12km s _1 , we derive = 220 ±11 km s _1 . Both values agree 
within their uncertainties, and averaging the two values we ob- 
tain = 214 + 7 km s _1 . With these values for R and O 
we can derive the systemic velocity v sys (local standard of rest 
velocity) measured from the Sun in the direction of Wd 1 as a 
function of distance (Fig. [TJ. 

We tested our assumption of a flat rotation curve by compar- 
ing the measured Tangent point velocity from the SGPS data 
with the velocity extrapolated from our values for . The 
Tangent point at any line of sight through the inner Galaxy is 
the point at which we look tangential to the circular motion and 
hence find the largest line of sight velocity. For flat rotation we 
should observe a radial velocity of -139 km s at the Tangent 
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Table 1. The systemic velocity and distance to the spiral arms 
along the line of sight toward Wd 1 determined with the Spiral 
arm model in Fig [2] and the rotation curve in Fig.Q] 



Spiral Arm 


v JVI [kms '] 


dm [kpc] 


Sagittarius I 


-17 


1.5 


Scutum-Crux I 


-52 


3.5 


Norma I 


-89 


4.9 


Norma II 


-35 


11.5 


Scutum-Crux II 


+6 


15.0 


Sagittarius II 


+26 


18.4 



point in the direction of / = 339.6° (Fig. [TJ. This agrees rea- 
sonably well with the observed value of somewhere between 
-130 km s _1 and -135 km s" 1 (Fig.0. This is actually a quite 
remarkable result, since the Tangent point in the direction of 
Wd 1 is at a distance of about 7 kpc from the Sun at a Galacto- 
centric radius of about 2.6 kpc. This should be well within the 
area which is heavily influe nced by the central bar of the Milky 
Way. Benjami n et alj d2005l) determined that the central bar has 
a maximum distance of 4.4 kpc from the Galactic centre, using 
a Galacto-centric distance of 8.5 kpc for the Sun. With the value 
of R G = 7.6 kpc this would translate to a maximum distance of 
3.9 kpc which is significantly higher than the Galactoc-centric 
distance of the Tangent point. 

To determine the distance to the spiral arms in the direction 
of Wd 1 we utilise the distribution of free electrons i n our G alaxy 
as described in the models of Tayl or & Cordel (1 19931) and 
ICordes & Laziol (|2002|) . The peaks in electron density indicate 
the locations of the spiral arms in this d irecti on of the Galaxy 
(Fig, pj i. Since iTavlor & Cordesl d 19931) and ICordes & Laziol 
d2002l) used a Galacto-centric radius of R G = 8.5 kpc for the 
Sun in their spiral arm model, we re-scaled this profile for 
R Q = 7.6 kpc. The distances to the nearby Hn regions used to 
define the distance of the nearby spiral arms were derived spec- 
troscopically. Hence, we left the distance to the two nearby spi- 
ral arms alone and only re-scaled the distances to the other more 
distant spiral arms. With the rotation curve in Figure [TJ we de- 
rive the systemic velocities at the centres of the individual spiral 
arms. These newly determined distances and systemic velocities 
are listed in Table [TJ A comparison of these values with the Hi 
emission profile in the direction of Wd 1 is shown in Fig. [3] 

The locations of the individual spiral arms in velocity space 
agree remarkably well with the brightness temperature peaks in 
the Hi emission profile (Fig. [3). The emission peak between 
-140 and -100 kms~' near the Tangent point is produced by a 
pile up in velocity space rather than a real density enhancement. 
In this region, the velocity does not change significantly with dis- 
tance, resulting in one velocity channel containing the emission 
from H i distributed over a much larger distance interval than in 
other locations along the line of sight. This creates a peak in the 
H i emission profile close to the Tangent point velocity. There 
are two spiral arms that can be seen distinct from the others, 
Norma I with a centre velocity of -89 km s _1 and Sagittarius II, 
for which the model centre velocity of +26 kms~' is very close 
to the emission peak at +22 kms -1 . Since Sagittarius II is sup- 
posed to be well outside the Solar Circle on the other side of the 
Galaxy in any Galactic model, this agreement with the H i emis- 
sion profile is incredible and boosts our confidence in our rota- 
tion model. The other spiral arms are clumped into two groups: 
Scutum-Crux I and Norma II merge at a radial velocity of about 
-40 kms -1 and the Local arm, Sagittarius I, and Scutum-Crux 
II are at low negative velocities. 
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Fig. 3. Hi emission profile averaged over an area of a few ar- 
cminutes towards Wd 1. The locations of the spiral arms as listed 
in TableQJare indicated. The Local arm is indicated at km s" 1 
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Fig. 4. H i absorption profiles of four radio bright H n regions in 
the vicinity of Wd 1 . To make the individual components bet- 
ter comparable we plotted the relative absorption - the ratio of 
absorbed emission to total emission - as a function of radial ve- 
locity. The dotted v ertical lines i ndicate the systemic velocities 
of the H ii regions (lRusseilll2003l) . which, within their uncertain- 
ties, are in excellent agreement with the absorption peaks. 



3.2. H i absorption 

The radio continuum emission of the H n region produced by 
the Wd 1 cluster is not bright enough at 1420 MHz to produce 
a prominent Hi absorption profile. Its peak brightness tempera- 
ture is only about 24 K in the SGPS data and the source is barely 
resolved. Its diameter is about 3' at a resolution of 2'. Therefore, 
we have to average over several velocity channels to amplify 
the weak absorption signal and increase its signal-to-noise ratio. 
For each spiral arm along the line of sight we tried to average 
those channels that are likely to produce the deepest absorption. 
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Table 2. Characteristics of radio bright H n regions and super- 
nova remnants in the vicinity of Wd 1. Column 1 contains a 
name based on the Galactic coordinates of the object. Columns 
2 and 3 contain information about the systemic velocity and de - 
rived kinematic distance found in the literature dRusseill 12003). 
Column 4 contains the distance to the source based on the rota- 
tion curve and method described in Sections 3.1 and 4. 1, respec- 
tively. Column 5 indicates the proposed spiral arm the object 
resides in. 



Source 


Vsys 


da, 


dnew 


Spiral Arm 




[kms- 1 ] 


[kpc] 


[kpc] 




G337.8-0.1 






all 


Norma II 


G337.95-0.48 


-41 


3.1 


2 9 +L2 


Scutum Crux I 


G338.0-0.1 


-51 


12.0 


10.7 +u 
14 + ^ 


Norma II 


G338.41-0.24 


-1 


15.7 


Scutum-Crux II 


G338.4+0.1 


-32 


13.1 


1 1 7 +20 

-0J 


Norma II 


G338.8+0.6 


-62 


4.3 


, q+oVsr 


Scutum-Crux I 


G338.5+0.1 






all 


Norma II 


G338.9-0.1 


-38 


3.1 


2-8^ 


Scutum-Crux I 


G338.9+0.4 






-3.9 


Scutum-Crux I 


G339. 13-0.41 


-38 


3.1 


36=3 


Scutum-Crux I 


Wd 1 


-55 




- ,4 °-n.4 


Scutum-Crux I 


G339.58-0.12 


-34 


2.8 


2 6 +13 

-on 


Scutum-Crux I 


G339.84+0.27 


-20 


14.1 


12 5 +2J 


Norma II 


G340.2-0.2 


-50 


3.7 


3.5+10 


Scutum-Crux I 


G340.24-0.48 


-61 


4.4 


q +8:§ 


Scutum-Crux I 


G340.6+0.3 






al5 


Scutum-Crux II 



In these absorption profiles we can identify those velocity in- 
tervals in the spiral arms where the H i produces deep and dis- 
tinct absorption features and those that do not. As a reference 
we used the absorption profiles of the nearby H n regions shown 
in Fig. E] Two of these, representing the G338.9+0.6 and the 
G340. 2-0.4 complexes, are Scutum-Crux I objects according to 
their systemic velocity (see Table |2). The other two Hn regions 
are located beyond the Tangent point in Norma II (G338. 4+0.1 
complex) and G338.41-0.24 might even reside in Scutum-Crux 
II. 

The absorption of the three nearby arms, the Local arm, 
Sagittarius I, and Scutum-Crux I, was identified in the absorption 
profiles of the two H n regions in Scutum-Crux I since their pro- 
files are not contaminated by absorption of H i located beyond 
the Tangent point. For Norma I and the Tangent point we used 
the two H ii regions that are located beyond the Tangent point 
on the other side of the Galaxy. We found distinct absorption in 
the Local arm between +5 and -1 kms -1 , for Sagittarius I be- 
tween -10 and -20kms _1 , for Scutum-Crux I between -32 and 
-45 km s _I , for Norma I between -85 and - 102 km s , and for 
the Tangent point between -110 and -125kms~'. For Scutum- 
Crux I we did not use negative velocities higher than -45 km s _1 
to avoid confusion with a possible bubble centered at Wd 1 (see 
section 3.3). 

The results of our averaging procedure are shown in Fig [5] It 
is apparent that the radio continuum emission from Wd 1 is ab- 
sorbed by the Local arm, Sagittarius I, and Scutum-Crux I. This 
is indicated by a hole in the distribution of H i emission at the po- 
sition of Wd 1 . There is no evidence for absorption in Norma I or 
the Tangent point gas. This indicates a location in Scutum-Crux 
I or at the near edge of Norma I. The latter, however, is rather 




339.6° 339.5° 

GALACTIC LONGITUDE 



Fig. 5. H i channels images averaged together to demonstrate in 
which spiral arms Wd 1 is absorbed. For the individual arms we 
averaged over: +5 to -1 kms _I (Local arm), -10 to -20 kms -1 
(Sagittarius arm), -32 to -45 km s _1 (Scutum-Crux arm), -85 to 
-102 kms -1 (Norma arm), and -110 to -125 kms -1 (Tangent 
Point). In these images white denotes high and black weak emis- 
sion. The radio continuum emission of Wd 1 is indicated by the 
white contours. 



unlikely, because there must not be any absorbing Norma I ma- 
terial between us and Wd 1 . All H n regions shown in Fig. |4]have 
a very deep distinct absorption feature at their systemic velocity, 
which is likely created by material in their vicinity in the cloud 
complex from which the stars in those H n regions were formed. 
Since the remains of these clouds are expected to be very dense 
and cold the optical depth should be rather high. This should pro- 
duce deep absorption lines. If Wd 1 is located in the Norma arm 
it cannot have such a component. This makes a Scutum-Crux 
arm location more likely, however, it does not entirely exclude 
the possibility of a location at the near edge of the Norma arm. 

3.3. The neutral environment of Wd 1 

We carefully investigated the HI data in the velocity range of 
the Scutum-Crux arm and at velocities corresponding to the near 
edge of the Norma arm. We found one prominent feature in ve- 
locity space which is exactly centered at the position of Wd 1 
at a velocity of -55 kms . This feature is shown in Fig [6] indi- 
cated by the dashed ellipses (top and right panels). The dashed 
circle in the centre image marks the location of this expanding 
bubble projected back to the map plane. This feature, Bl, cannot 
be an absorption feature since the depression visible at the loca- 
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I = 339.56° 




v = —55 km s 



339.7° 339.6° 339.5° 339.4° 
Galactic Logitude 

Fig. 6. Three views of the two small expanding bubbles, Bl and B2, around Wd 1, indicated by dashed and dotted ellipses respec- 
tively. Longitude-velocity slices are displayed at the top and bottom and the velocity-latitude slices in the left and right panels. The 
HI channel map at about -55 km/s in longitude-latitude is shown in the centre, with Wd 1 indicated by the white contours. The 
bright emission in the velocity slices mark the gas at the far edge of Scutum-Crux I. In these images black represents low and white 
high intensity. 



tion of Wd 1 is deeper than the peak brightness of Wd l's radio 
continuum emission. In the absorption map of the Scutum-Crux 
arm (Pig. [5j the velocity interval related to this expanding bub- 
ble was omitted to avoid confusion with a hole in the H i map 
that is created by a lack of H i and not actual H i absorption. 

Feature B 1 seems to be located at the edge of a much larger 
elliptical bubble (B2), indicated by the dotted ellipses in Fig. [6] 
This feature is very obvious in the map plane but is not as clear 
in velocity space (left and bottom panels of Fig. |5). There are 
fingers of emission emerging from the bright Scutum-Crux I gas 
towards higher negative velocities, indicating an expansion for 
this feature. However, these fingers are not closed in an end cap, 
which would mark the part of the bubble that is moving towards 
us. Bubble Bl is closed in velocity space, which makes it very 
easy to detect. The expansion velocity of B 1 is 5 km s~ 1 , as read- 
ily seen in Fig. [6] The expansion velocity of B2 is not easy to 
identify since the end cap is missing, but the velocity fingers 
seem to indicate a somewhat higher expansion velocity. The spa- 
tial coincidence of the two features suggests both are associated 
with Wd 1, supported by their central velocity of -55 kms -1 . 
This velocity is well within the velocity interval predicted for 
Wd 1 by the Hi absorption measurements. Therefore we adopt 
a radial velocity of -55 + 3 kms' 1 for Wd 1. Uncertainties for 
the radial velocities along the line of sight are discussed in Sect. 
4.1. 

At a radial velocity of -55 kms 1 a larger field of view re- 
veals a much larger bubble-like feature (B3), which is open to 
the south, away from the Galactic plane (Fig. |7J. The emission 



associated with Bubble B2 can be seen in the lower right corner 
of B3. To the north, B3 consists of a shell with two large, bright, 
and complex emission regions to the east and west. The eastern 
region contains the Hn complex G340.2-0.2, with a radial ve- 
locity of -50 kms -1 (Tabled. This suggests that Wd 1 and the 
Hn region complex G340.2-0.2 are evolving in the same envi- 
ronment, at the far side of Scutum-Crux I. 



4. Discussion 

4. 1 . The distance to Westerlund 1 

We determine a radial velocity of -55 ± 3 km s _1 for Wd 1 
from the central velocity of the bubbles in which we believe 
it is located. Using the rotation curve determined in Sect. 3.1 
(see Fig. [TJ this translates to a distance of 3.6 + 0.2 kpc or 
10.6 + 0.2 kpc. The closer distance is preferable, because the H n 
region of Wd 1 does not indicate any absorption in the Hi gas 
of Norma I or near the Tangent point. Additionally, the previous 
distance estimates are ~ 4 - 5 kpc. A comparison with the spi- 
ral arm model in Fig.|2]and the calculated distances to the spiral 
arms in Table Q~]reveals that Wd 1 is located in the Scutum-Crux 
I arm. 

A possible source of uncertainty for our rotation curve could 
be the presence of non-circular motion. Non-circular velocity 
components due t o stream ing motion can reach a maximum of 
about +10 km s _1 dBurtonl 19721) . This would result in +8 km s _1 
projected to the line of sight towards Wd 1 assuming a distance 
of 3.6 kpc. Another source of non-circular motion could be a 
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340.0° 339.5° 

GALACTIC LONGITUDE 



Fig. 7. H i channelmap at a radial velocity of -55 km s 1 from the 
SGPS. Black contours go from 20 K to 90 K in steps of 10 K. 
Grayish contours indicate the 1420 MHz continuum emission. 
The location of Wd 1 is indicated. 



spiral shock, which is always directed towards the Galactic cen- 
tre. At the location of Wd 1 the spiral shock would be directed 
away from us, leading to an underestimate of the actual distance. 
Typical values for the velocity shi ft due to the spiral shock coul d 
reach a maximum of 30 km s^ 1 dFoster & MacWilliamsi r2006). 
which would translate to about 24 km s projected to the line of 
sight. Random motion add another ±5 km s to our uncertainty. 
Taking all these uncertainties into account, we derive a systemic 
velocity of km s _1 forWd 1. In the direction of Wd l,the 

velocity versus distance gradient is rather steep (see Fig [TJ, so 
the large velocity uncertainty does not translate to a large dis- 
tance uncertainty. This radial velocity corresponds to a distance 
of 3.6*Q4 kpc, or by choosing the centre of the uncertainty in- 
terval as the most probable distance, we get 3.9 ± 0.7 kpc. This 
result agrees very well with the latest indepe ndent distance esti- 
mates of 4.0 + 0.3 kpc dBrandner et al.ll2005l) and 4.7 ± 1.1 kpc 
dCrowther et al.l 120061) . This distance is consistent with a loca- 
tion at the far side of the Scutum-Crux I arm. Furthermore, the 
bubbles Bl and B2 in Fig.[6]seem to be emerging from the bright 
gas at the high velocity edge of the Scutum-Crux I arm towards 
higher negative velocities, supporting a location on the far side 
of this arm. 



4.2. Distance estimates for H n regions and SNRs in the 
vicinity of Wd 1 

In Figure [8] we show a 1420-MHz radio continuum image of 
the area around Wd 1, created from the end-channels of the Hi 
data that do not contain any neutral hydrogen emission. Wd 1 
is located in a rich area of radio bright H n regions and super- 
nova remnants. In Table [2] we list H n regions and complexes for 
which the systemic velocity is known and for which the distance 
ambiguity in the inner Galaxy, where each radi al velocity cor - 
responds to two possible distances, was solved dRusseilll2003l) . 
With the rotation curve shown in Fig.Q] new distances to these 
objects can be deduced. A comparison of the systemic velocities 
of these H n regions with the centre velocity of each spiral arm, 
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Fig. 8. Radio continuum emission at 1420 MHz around Wd 1 
taken from the end channels of the SGPS Hi observations. 
Contours are at 11, 17, 23, and 40K (black) and 80, 120, and 
160 K (white). The dotted ellipses enclose SNRs and H n regions 
that are believed to be beyond the Tangent point. The location of 
Wd 1 is indicated. 
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Hi absorption profiles of one radio bright Hn region 
9+0.4) and three bright SNRs in the vicinity of Wd 1. 



determined in Sect. 3.1 (see Table [T), gives the most likely spi- 
ral arm where these objects reside (Table|2]i. Almost all of these 
sources are apparently concentrated in the Scutum-Crux I arm 
and the Norma II arm. Only G338.41 - 0.24 seems to be located 
even further away in the Scutum-Crux II arm. Objects believed 
to be beyond the Tangent point are encircled by dashed ellipses 
in Fig. [8] It is quite curious that no objects were found in either 
Sagittarius I or Norma I arms 

Hi absorption profiles for four additional radio bright 
sources can also be determined (Fig. |9). The distances and 
systemic velocities of these objects were previously unknown. 
These absorption profiles were compared with those in Fig|4]to 
determine the most likely spiral arm location for these objects. 
We then assume the distance to the centre of that particular spiral 
arm to be the most likely value for these objects, with uncertain- 
ties ~ ±1 kpc (Table 0. 

The absorption profile of the HII region G338.9+0.4 in Fig.|9] 
is very similar to G338.9 + 0.6 and G340.2 - 0.2 (see Fig.gJ with 
no evidence of absorption in Norma I or at the Tangent point. 
This suggests G338.9 + 0.4 could be a part of the G338.9 + 0.6 
H n region complex, just to the north (see Fig [8]). Therefore, we 
propose a Scutum-Crux I location at a distance of about 3.9 kpc. 
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Fig. 10. Radio continuum spectrum of the H n region, produced 
by the members of Wd 1 . 



The SNR G340.6+0.3 is absorbed by Norma I and the 
Tangent point gas and also shows a deep absorption feature with 
a relative absorption of almost 0.7 at about km s _1 (Fig[9j». This 
makes its absorption profile very similar to that of G338.41-0.24 
(Fig |4j. No other source shows such a deep and wide absorp- 
tion feature around a velocity of kms . This indicates a loca- 
tion on the other side of the Galaxy in Scutum-Crux II at about 
15 kpc. This gives the SNR a diameter of about 26 pc. It is not 
unusual that absorption features that are far away seem to be 
deeper than those that are nearby. The reason for this is simply 
the area probed by the observing beam becomes bigger with dis- 
tance so that more H i is detected within one beam. 

The SNRs G337.8-0.1 and G338.5+0.1 show absorption by 
Norma I and at the Tangent point but lack the deep and wide 
absorption feature at km s _1 . Their absorption profiles look 
very similar to that of G338.4+0. 1 in Fig. [4] Therefore these two 
SNRs are most likely Norma II objects at a distance of about 
11 kpc. This gives G338. 5+0.1 a diameter of about 30 pc and 
G337 . 8 -0. 1 an extent of about 30x20 pc. 

4.3. The Wd 1 H n region 

We determined radio continuum flux densities for the Wd 1 H n 
region at 1 .4 GHz (620 ± 60 mJy) from the continuum end chan- 
nels of the SGPS Hi data and a t 843 MHz (680 + 80 mJy) 
from data of the Molonglo Survey dGreen et al.f l999). We added 
two flux de nsity measurements from the literature at 5 GHz 



(520 mJy, 



Wright et al. 



Havnes et alJ fl979l) and at 4.85 GHz (474 mJy, 



1994) (Fig. [10>. Since no error was given for the 



Havnes et al. I ( 1 19791) value we weighted both fluxes equally. 



The error bars plotted in Fig. [TO] represent 10% of the flux 
density. A weighted least-squares fit gives a spectral index of 
a = -0.17 + 0.07, that we interpret as arising from optically thin 
thermal plasma. 

The observed diameter of Wd 1 in our 1.4 GHz data is 3'.5 at 
a resolution of 2'. This results in an actual diameter of 3'. This 
agrees with its app earance in the 843 MHz Molonglo Survey 
dGreen et al.| [l999). The observed peak of the radio continuum 
emission is about 24 K. Correcting for the synthesised beam, the 
H ii region has a peak brightness temperature of 33 K at 1 .4 GHz. 
The observed brightness temperature Tb at frequency v is related 



to the electron temperature T e inside the H n region and the op- 
tical depth r by: 



T b (v) = T e (l-e-™), 

where the optical depth at frequency v is given by 
'TeV 1 -* 3 I v r 21 / EM 



*> - 8.235 



pc cm 



-6 



(1) 



(2) 



(e.g. lRohlfs & Wilsonll2004l) . Here EM is the emission measure 
related to the path length / through the H n region and the elec- 
tron density n e inside it by: 



EM 



dl 



(3) 



With an electron temperature of T e - 8000 ± 2000 K and a path 
length of 3.5 pc (3' at 3.9 kpc) through the Hn region, we esti- 
mate an optical depth of 0.0041 ^J"^ and an emission measure 



of 19000!*°™ pc cm- 6 



Hence, n e — 74 ± 13cirT 3 and the total 
mass of the ionized gas inside the Hn region is 53 ± 9M . The 
frequency at which this emission becomes optically thick (i.e. 
t v = 1) is ~ 100 MHz. 

Alternatively, the electron density of the H n region can be 
calculated from the observed flux density. For an optically thin 
Maxwellian plasma, the flux S v is given by: 



S v = 5.7 10 -56 T;°- 5 g ff E v d- 2 mJy, 
where gfj is the Gaunt-factor defined by: 



Sff = 



V3 



17.7 + In 



and Ey is the volume emissivity: 
E v = j n 2 e dV. 



(4) 



(5) 



(6) 



The radio spectrum in Fig. [10] indicates that the emission is 
optically thin over the entire observed frequency range. Using 
the flux density at 843 MHz, we derive an electron density of 
n e = 65 + 7 cirT 3 and a total ionized mass of M — 47 ± 5 M Q . 
This is in excellent agreement with the value derived from the 
observed brightness temperature. 

4.4. Wd 1 and its neutral environment 

In Section 3, two small (Bl and B2) and one large (B3) bub- 
ble in the H i data were described. We believe these are related 
to Wd 1. The determination of the mass of atomic material in 
these bubbles is difficult since these objects, in particular bubble 
B3, seems to merge with the surrounding ISM in the Hi images 
( Figs. |6] and |7j. Bubble Bl and B2 together contain a mass of 
about 300 M with an error of at least 50 %, assuming the H i is 
optically thin and the mass ratio between hydrogen and helium 
is 10:3. The inner edge of bubble Bl seems to be just outside the 
Wd 1 Hn region (see Fig.|6]l giving a diameter of about 5 pc at a 
distance of 3.9 kpc. At t he same distance, bubble B 2 h as an ex- 
tent of about 18x10 pc. iMcClure-Griffiths et ail d2002l) defined 
the dynamic age of a stellar wind bubble that expands as oc f° 3 
by: ' 



t 6 = 0.29 , 



(7) 
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where ?6 is the dynamic age of the bubble in Myr, R the radius 
in pc, and v exp the expansion velocity in kms -1 . For bubble B2 
the maximum expansion radius calculated from the position of 
Wd 1 is about 10 pc. If we use an expansion velocity of 5 km s _1 
for B2, determined for bubble Bl, we derive a dynamic age of 
about 600,000 yr for B2. Since its structure in velocity space (see 
Fig. [6]i implies a somewhat higher expansion velocity than that 
for bubble Bl and the radius is a maximum expansion radius, 
this a ge represents an upp er limit. 

ICrowther et all ([2006) determined an age of 4.5 - 5.0 Myrs 
for Wd 1, based on the ratio of WR stars to red and yellow hy- 
pergiants. Clearly, bubbles Bl and B2 cannot be the stellar wind 
bubble produced in the early evolution of Wd 1 ; they are sim- 
ply too young. Both Bl and B2 could have been created after 
the last supernova explosion pushed away all material inside the 
stellar wind bubble. The recent discovery of the X-ray pulsar 
CXO J 1647 10.2-4552 16 inside Wd 1 imp lies that at least one 
supernova explosion occured inside Wd 1 (Muno et al. 2006|). In 
this case we speculate that Bl and B2 originate from the only 
source of matter and energy left after the last supernova oc- 
curred, the combined winds and mass loss of the cluster stars. 
In that scenario, the gas in the shells of those bubbles must con- 
sist mostly of recombined material previously ejected by the 
stars in their winds, since the supernova shock wave should have 
removed any remaining material in the ambient medium. This 
would imply that the last supernova explosion happened less 
than 600,000 yr ago. Additionally, the ionized material inside 
the H ii region would also be material ejected by the stars in their 
winds. The total mass of the ionized gas plus the B2 Hi shell is 
~ 350 M Q . If we assume these structures are 500,000 yr old we 
require 70 stars with an average mass-loss rate of 10~ 5 M Q yr 
for each star, which is highly feasible for the stellar population 
inWdl. 

We suspect that the impact of the original stellar wind bub- 
ble, created by the members of Wd 1 in their early stage of evo- 
lution, is represented by bubble B3 (Fig. [7]). Its diameter is about 
50 pc at a distance of 3.9 kpc. The location of Wd 1 in the lower 
right corner of this feature is likely the result of a highly struc- 
tured environment. To the South, away from the Galactic plane, 
the density is expected to decrease in which case the winds of the 
members of Wd 1 pushed the material away from the Galactic 
plane into the Galactic halo. To the West we find a dense cloud 
complex that impeded free expansion, whereas to the East the 
winds of the stars pushed material into the cloud centered at 
the G340.2 - 0.2 Hn region complex. In Fig. [8] we can identify 
three bright compact radio sources that belong to this complex 
at positions (Lb): (340.3,-0.2), (340.1,-0.15), and (340.1,-0.25). 
These are all sources with far infrared colours characteristic of 
ultra-compact Hn r egions, which are kno wn to be regions of 
OB star formation (Bronfm an et al.| [l996). This indicates that 
the G340.2-0.2 Hn region complex is younger than Wd 1 and 
its formation may well have been triggered by the interaction of 
the early Wd 1 stellar wind or a supernova explosion with the 
clouds in which G340.2-0.2 is embedded. 

To the North of bubble B3 we find a shell that might still 
be expanding. Since this shell is the only dynamic feature we 
observe, we use its distance of 45 pc from Wd 1 as the expansion 
radius. We do not see any expansion of B3 in the velocity slices 
of our H i data. If we assume that the dynamic age of bubble B3 
equals the assumed age of Wd 1, an expansion velocity of about 
3 kms -1 is implied, which is a reasonable value. Bubble B3 is 
certainly much older than B2 and Bl, because it is many times 
as large and a low velocity is indicated by the absence of any 
visible feature in the velocity slices. Unfortunately we cannot 



estimate the mass of the bubble because of the confusion with 
the emission regions to the east and the west. 

5. Summary 

Based on Hi observations from the SGPS, we have established 
that Wd 1 is in the Scutum-Crux arm of the Galaxy. We find 
three large bubbles in which Wd 1 is located at a radial velocity 
of -55^ 6 km s _1 . Using a flat rotation model of the Galaxy and 
adopt ing a galacto-centric distance of 7.6 kpc dEisenhauer et alj 
2005) and a velocity of 214 + 7 kms' 1 for the Sun around the 
Galactic Centre, we transpose the radial velocity of the bubble 
features to a distance of 3.9 ± 0.7 kpc at a Galacto-centric ra- 
dius of about 4.2 kpc. We are confident in our derived distance- 
velocity calibration in this direction, since it predicts to within 
+4 kms~' the velocity of the Tangent point and the velocity of 
the H i gas in the Sagittarius Arm outside the Solar Circle on the 
far side of the Galactic Centre. Moreover, this is very encourag- 
ing since it is believed that the position of Wd 1 is in the region 
of the Galactic bar. 

A distance of 3.9 ± 0.7 kpc is somewhat less than the ear- 
lier estimate of 4.7 ± 1. 1 kpc based on the a bsolute magnitude 
calibration of WR stars ( Crowther et al. 2006), but very s imilar 
to the 4.0 + 0.3 kpc determined by iBrandner et alj (|2005) from 
initial analysis of the photometric detection of the main and pre- 
main sequence populations in Wd 1 . Though the uncertainty on 
these estimates is relatively large (~ 10 — 25%), a weighted aver- 
age of these values gives us confidence that the distance to Wd 1 
is ~ 4 kpc. Until t he final analysis of t he photometric observa- 
tions described in Bran dner et al.l d2005l) is available, we argue 
that this represents the best distance estimate to Wd 1 . 

A study of the H i images from the vicinity of Wd 1 revealed 
three bubble-like features Bl, B2, and B3. We believe that bub- 
bles Bl and B2 are the result of the stellar winds of the members 
of Wd 1 and contain recombined wind material. We argue bubble 
B3 represents the stellar wind bubble created in the early history 
of Wd 1 . The formation of the H n region complex G340. 2-0.2 
which is embedded in the eastern edge of B3 may have been trig- 
gered by the interaction of the Wd 1 wind with the dense clouds 
in this area. If this is true, the Wd 1 G340.2-0.2 pair would be a 
nice example of sequential star formation. 
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